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Abstract 

Since Grover's seminal work which provides a way to speed up combinatorial search, 
quantum search has been studied in great detail. We propose a new method for designing 
quantum search algorithms for finding a marked element in the state space of a graph. 
The algorithm is based on a local adiabatic evolution of the Hamiltonian associated with 
the graph. The main new idea is to apply some techniques such as Krylov subspace 
projection methods, Lanczos algorithm and spectral distribution methods. Indeed, using 
these techniques together with the second-order perturbation theory, we give a system- 
atic method for calculating the approximate search time at which the marked state can 
be reached. That is, for any undirected regular connected graph which is considered 
as the state space of the database, the introduced algorithm provides a systematic and 
programmable way for evaluation of the search time, in terms of the corresponding graph 
polynomials. 

Keywords: Quantum search algorithm, Local adiabatic evolution, Graph, 
Krylov subspace, Lanczos algorithm, Spectral distribution, Second-order per- 
turbation theory. 
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Grover's quantum search algorithm [1] is one of the main applications of quantum computation. 
This algorithm is sometimes described as a way for searching a marked item in an unstructured 
database of N items in time 0(y/~N). This gives a quadratic speedup over the exhaustive 
search for a variety of search problems [2]. But the algorithm as originally proposed is not 
designed to search a physical database. In [3], Aaronson and Ambainis present a model of 
query complexity on graphs, where they showed that a database of iV items laid out in d 
spatial dimensions can be searched in time of order y/~N. In [4], Childs and Goldstone have 
considered an alternative quantum search algorithm based on a continuous time quantum 
walk on a graph. Quantum walks provide a natural framework for the spatial search problem 
because the graph can be used to model the locality of the database. In fact, for the case 
of the complete graph (unsorted database), the resulting algorithm is simply the continuous 
time search algorithm of Farhi and Gutmann [5]. On the hypercube, their results showed that 
the algorithm also provides quadratic speedup [6, 7, 8]. In [6], Farhi et al. have used the 
time-dependant Hamiltonian approach for Grover's problem, where they have considered the 
constant-rate adiabatic approach (the requirement of adiabaticity is expressed globally) and 
obtained a complexity of order N, the number of items in database. Then, Roland and Cerf 
[9] have been considered the same problem with the same approach with this deference that 
they considered the adiabatic evolution locally, i.e., they adjusted the evolution rate of the 
Hamiltonian so as to keep the evolution adiabatic on each infinitesimal time interval and by 
this adjusting the total running time of order \/N has obtained. Recently, another adiabatic 
version of the quantum search problem has considered by Rezakhani et al. in [10], where they 
have employed continuous time, global adiabatic evolution in order to calculate a quantity 
called "adiabatic error" , which quantifies the distance between the instantaneous ground state 
and the actual marked state. 
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In this paper, we follow the approach of the paper [9] and employ some techniques such 
as Krylov subspace projection methods, Lanczos algorithm, spectral distribution methods and 
the second-order perturbation theory in order to give a systematic method for calculating the 
total search time approximately. In fact, by reducing the Hilbert space of the corresponding 
Hamiltonian to the smaller one called Krylov subspace and using the spectral methods, we 
calculate the minimum energy gap between two lowest eigenvalues of the Hamiltonian and 
consequently the total search time, in terms of the polynomials associated with the graph. The 
organization of the paper is as follows: In section 2, we recall some preliminary facts about 
Krylov subspace projection methods and spectral distribution method needed for the approach 
of the paper. In section 3, by using the local adiabatic evolution approach, quantum search on 
a graph is investigated, where an analytical but approximate formula for the quantum search 
time (the main result of the paper) is obtained in terms of the graph polynomials. Section 4 
is devoted to some examples of graphs, in order to clarify the introduced method in details. 
The last section contains a brief conclusion. 

2 Krylov subspace and spectral methods 

In this section we give a brief review of some of the main features of Krylov subspace projection 
methods and spectral distribution method (the reader is referred to Refs. [11]- [19] for more 
details). 

2.1 Krylov subspace projection methods 

Krylov subspace projection methods (KSPM) are probably the most important class of projec- 
tion methods for linear systems and for eigenvalue problems. In KSPM, approximations to the 
desired eigenpairs of an n x n matrix A are extracted from a d- dimensional Krylov subspace 
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where |0o) is often a randomly chosen starting vector called reference state and d <^ n. In 
practice, the retrieval of desired spectral information is accomplished by constructing an or- 
thonormal basis V d G R nxd of K d (\<f> ),A) and computing eigenvalues and eigenvectors of the 
d by d projected matrix H d = Pv d T APv d , where Py A is projection operator to <i-dimensional 
subspace spanned by the basis V d . 

The most popular algorithm for finding an orthonormal basis for the Krylov subspace, is 
Lanczos algorithm. The Lanczos algorithm transforms a Hermitian matrix A into a tridiagonal 
form iteratively, i.e., the matrix A will be of tridiagonal form in the <i-dimensional projected 
subspace H d . In fact, the Lanczos algorithm is deeply rooted in the theory of orthogonal 
polynomials, which builds an orthonormal sequence of vectors (|</>o), |</>i), • \4>d-i)} and satisfy 
the following three-term recursion relations 



A\<f>i) = /3 i+1 \(j) i+1 ) + Oilfc) + 



(2-2) 



The vectors \4>i),i = 0,1,. ..,d — 1 form an orthonormal basis for the Krylov subspace 
K d (\(p ),A). In these basis, the matrix A is projected to the following symmetric tridiagonal 
matrix: 

/ 

a ft 
h a, (3 2 
(3 2 a 2 (3 3 

(2-3) 



« d -2 Pd-i 



V 



&_ 2 

p d -i Oi d -i J 

where the scalars f3 i+ \ and ctj are computed to satisfy two requirements, namely that \4>i+i) be 
orthogonal to and that |||0 i+ i)|| = 1. 

In fact, the Lanczos algorithm is a modified version of the classical Gram-Schmidt orthog- 
onalization process. As it can be seen, at its heart is an efficient three-term recursion relation 
which arises because the matrix A is real and symmetric. 
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If we define the Krylov matrix K such that the columns of K are Krylov basis {A l (f)o; i = 
0, d — 1} as: 

ir:=(|0o),^|0o),...,^- 1 |0o)), 
the application of the orthonormalization process to the Krylov matrix is equivalent to the 
construction of an upper triangular matrix P such that the resulting sequence <3> = KP satisfies 
$t$ = 1. We denote by |0j) and Pj respectively the j-th column of $ and P. Then we have 

(001^(^)^(^)100) = {KP t \KP 3 ) = (0,10,), (2-4) 

where Pi = ao + a\A + ... + aiA 1 is a polynomial of degree % in indeterminate A. 

In the remaining part of this section we give an algorithmic outline of the Lanczos algorithm, 
where it will be used in calculation of parameters a, and appeared in (2-2). 
Lanczos algorithm 

Input: Matrix A G R nxn , starting vector |0 O ), |||0o)|| — 1, scalar d 
Output: Orthogonal basis {|0o), |0d-i)} °f Krylov subspace K d (\(f) ), A) 

A) = o,|0o) = |0)/|||0)|| 

for 2 = 0,1,2,... 

\vi) = A\(j) t } 

Oti = (<f>i\Vi) 
\v i+1 ) = \Vi) - &|0i_i) - Oi\<j>i) 

A+i = llk+i)ll 

if 

A+i ^ o 

|0 i+ i) = 

= o. 
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2.2 Spectral distribution method 

For any pair (A, |0 O )) of a matrix A and a vector \<f>o), one can assign a measure \i as follows 

Ai(a;) = <0o|£(a;)|0o), (2-5) 

where -E(x) = |wj)(tii| is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE(x)dx (2-6) 
so that, for any polynomial P(A) we have 

P(A) = J P(x)E(x)dx, (2-7) 

where for discrete spectrum the above integrals are replaced by summation. For example, 
in continuous time quantum walk on a graph [21], the expectation value of powers of the 
corresponding adjacency matrix A over starting site |0o) can be written as 

(0o|^ m |0o) = / x m fi(dx), m = 0,1,2,.... (2-8) 

JR 

The existence of a spectral distribution satisfying (2-8) is a consequence of Hamburgers theo- 
rem, see e.g., Shohat and Tamarkin [[16], Theorem 1.2]. 

From the orthogonality of vectors \<f)j) (Krylov bases) produced from Lanczos algorithm 
process we have, 

*y = = (MPh^AMo) 

= J P*(x)P j (xMx)dx = (Pi.P.-V (2-9) 

Conversely if Po, Pd-i is the system of ortho normal polynomials for the measure \i then the 
vectors 

|0,) = P,(A)|0 o ), (2-10) 

will coincide with the sequence of orthonormal vectors produced by the Lanczos algorithm 
applied to (A, \<j>o)). 
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Now, substituting (2-10) in (2-2), we get three term recursion relations between polynomials 
Pj(A), which leads to the following three term recursion between polynomials Pj(x) 



(3 k +iPk+i(x) = {x- a k )P k {x) - P k Pk-i{x) 



(2-11) 



for k — 0, d — 1. Multiplying by f3\...(3 k we obtain 



f3 1 ...f3 k+1 P k+1 (x) = (x- a k )p 1 ...p k P k (x) - (3 2 k .(3 1 ...(3 k „ 1 P k - 1 (x). 



(2-12) 



By rescaling P k as Q k = (3i...f3 k P k , the three term recursion relations (2-11) are replaced by 

Qo(x) = l, Q 1 (x)=x, 
Qk+i{x) = (x - a k )Q k (x) - filQk^x), (2-13) 

for k — 1, 2, d. 

In the next section, we will need the distinct eigenvalues of adjacency matrix of a given 
undirected graph which can be written in the form (2-3), and the corresponding eigenvectors 
in order to obtain the minimum time at which the marked state can be reached. As it is known 
from spectral theory, we have the eigenvalues Xi of the adjacency matrix A as roots of the last 
polynomial Qd+i(x) in (2-13), and the normalized eigenvectors as [18, 19] 

/ Po(xi) X 



Pl(Xi) 



(2-14) 



3 Quantum search via local adiabatic evolution 

In this section, we investigate quantum search in an structured database by using the time- 
dependent Hamiltonian approach and spectral distribution method. To do so, first, we recall 
briefly the time-dependent Hamiltonian approach and local adiabatic evolution employed in 
Ref.[9] to Grover's problem. 
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Consider the evolution of a quantum system state \ip(t)) subject to a time-dependent Hamil- 
tonian H(t) is described by the Schrodinger equation (% = 1) 

ij t \m) = H(t)\m)- (3-15) 

According to the adiabatic theorem [20] , the state of the system will stay close to the instanta- 
neous ground state of the Hamiltonian at each time t, if the Hamiltonian varies slowly enough. 
In other words, for the instantaneous energy eigenbasis defined by H (t)\E n (t)) = E n (t)\E n (t)) , 
if we define the minimum gap between the lowest two eigenvalues as 

g mm = mi^iE^t) - E (t)] (3-16) 

and the maximum value of the matrix element of ^ between the two corresponding eigenstates 
as 

dH 

D_=max|(— > 1>0 | (3-17) 

with (^f-)i,o = {Ei(t)\^-\E (t)} , then the adiabatic evolution theorem guaranties that for the 
system which is prepared in its ground state \E (0)) at time t = 0, and evolved under the 
Hamiltonian H{t), we will have 

|(£ (TMT)>| 2 >l-e 2 (3-18) 

provided that 

< e, (3-19) 

9min 

where e C 1. 

3.1 Quantum search in an unsorted database with local adiabatic 
evolution 

Assuming a set of items in an unsorted database (uniform distribution) where one of them 
is marked, the main goal is finding the marked item in a minimum run time. To this end, a 
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quantum state \i) is assigned to each item i, so that the state space of the database is spanned 
by the states \i) with % — 0, 1, . . . , N — 1, where the marked state is denoted by \m). Since \m) 
is unknown a priori, the initial state is chosen as an equal superposition of all basis states, i.e., 

N-l 

Then, the Hamiltonian of the system is initially chosen as 



l^o> = 4= E l<>- (3-20) 



H = I N -\il)o){iPo\ (3-21) 

whose ground state is |-?/>o) with energy zero. Since the database is unsorted, we can associate 
the items of the database with the nodes of the complete graph with N vertices (denoted 
by Kn), for which the adjacency matrix is given by A = Jn — In- Therefore the initial 
Hamiltonian can be written as 

-. N-l ii i 

H o = In - „ E Q 1001 =In-j^Jn = ^[(N- 1)I n -A} = -L, (3-22) 

where, J N is an N x N all one matrix and L = (N — 1)In — A is the Laplacian of the complete 
graph Kn- Now, suppose that we are able to apply to our system the Hamiltonian 

H m = I N - \m)(m\, (3-23) 

whose ground state is the unknown marked state \m). Then, the corresponding time dependent 
Hamiltonian is considered as the following linear interpolation between H and H m 

H(s) = (1 - s)H + sH m = ^[(N-l + s)I N - (1 - s)A - Ns\m)(m\], (3-24) 

where s is an evolution function of time t and must be found optimally with the boundary 
conditions s(0) =0 and s(T) — 1. In the global adiabatic evolution s(t) is chosen as a linear 
function of t as s(t) = t/T, but in the local adiabatic evolution the time interval T is divided 
into infinitesimal time intervals dt and the adiabaticity condition is locally applied to each 
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of these intervals. Applying Eq.(3-19) to each infinitesimal time intervals, the local adiabatic 
condition is given by 



" S| <^ (3-25) 



dt 



for all times t. Now, by using the fact that |(^f-)i,o| < 1, the condition (3-25) is verified 
provided that 

\%\=^ 2 {s). (3-26) 

The algorithm consists in preparing the system in the state \ip ) and then applying the 
Hamiltonian H(s) during a time T so that s(T) = 1. In order to obtain the eigenvalues of 
H(s) and evaluate (3-16), we use the Lanczos algorithm to obtain the orthonormal basis in 
which the hamiltonian H(s) can be reduced to as a tridiagonal matrix. To do so, we choose 
\m) = |0 O ) as the reference state (starting vector in the Lanczos iteration algorithm) and apply 
H(s) on it to obtain the vector orthogonal to |0 O ) as 

l<M = — =(VjV|lfo> - |0o». 

Now, the Hamiltonian H(s) can be represented in the new orthonormal basis states |0 O ) an d 
|0x) as follows 

(l-s)(iV-l) -y/N^l(l - s) ^ 
-y/N=T(l-s) l + s(N-l) j 



\ 



(3-27) 



The eigenvalues of H(s) are given by E± = ^{N ± ^N 2 (l - 2s) 2 + 4iVs(l - s)}, so that the 
difference between these two eigenvalues gives the gap g(s) as 



N-4s(l -s)(N- 1) 

9(*) = y (3-28) 

We see that the minimum gap g min = ^= is attained for s = 1/2. Now, by using the local 
adiabatic condition (3-26) we obtain 

1^1 = ^-4^(1-^-1)]. (3-29) 



Structured quantum search 12 

After integration, one can find 

1 N 



t = — -j= = {arctanfViV - 1 (2s - 1)] + arctan y/N - 1}. (3-30) 

One may now evaluate the computation time of the algorithm by taking s — 1. For N 1, 
we obtain 

T ~ ^-VN, (3-31) 

which gives a quadratic speed-up with respect to a classical search, so that the algorithm can 
be viewed as the adiabatic evolution version of Grover's algorithm. 

3.2 Quantum search in a structured (non-uniform) database 

In the case that our states distributed non-uniformly, the state space of the database can 
be considered as an arbitrary connected graph other than the complete graph. Therefore, a 
straightforward generalization of the relations (3-22)-(3-24) leads us to consider 

H = 7 L, H m = I — \m)(m\ (3-32) 

so that 

H(s) = 7(1 - s)L + s(I - \m)(m\) = aI N - 7(1 - s)A - s\m)(m\, (3-33) 

where L is the Laplacian of the graph and 7 is a constant parameter which is determined in such 
a way that the search time be minimum (the search algorithm be optimal). We will consider 
regular undirected graphs so that the corresponding adjacency matrices are symmetric, and 
moreover the Laplacian L and the adjacency matrix A of the graphs differ from each other in 
only multiple of identity matrix, i.e., we have L = DI — A with D as the degree of each vertex. 
So, the parameter a in (3-33) is given by a = jD(l — s) + s. Since we need the minimum gap 
between two lowest eigenvalues of the hamiltonian H(s), so the first term in (3-33) (multiple 
of identity) can be dropped. 
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By choosing \m) = |0o) as the starting state for the Lanczos iteration algorithm, the 
corresponding Krylov bases are obtained via the three term recursion relation (2-2). That is 
we have, 

H(s)\(h) = (l(s - l)oo - s)\(fo) + 7 (s - 
ff(s)|0 i )= 7 (s-l){A|0i-i)+ai|0i)+A+i|0i+i)} ; i = l,2,...,d. (3-34) 
In other words, in the krylov bases \<f>i), the Hamiltonian H(s) is rewritten as 



(3-35) 



where, 



H. = A = 



a ft 

/3i ai & 

ft ••■ ••• 

o o ••. ■•. 








(2-1 



#1 



1 ... 



; o o 

: '•• 

... 



\ ... (3 d -! a d -i j 
Now, using (2-14) and the form of Hi, one can easily evaluate 

1 



so that we obtain 

= mHi\^) 



E 



(2) 



E 



k^i^i|^>i 2 _ i 



E 



i 



... 

... 

'•• 

'•• 





t+i Ef ] - Ef ~ l(s - 1) ( Xi - Xj ) Et v=0 IfMP&ixj) ' 
Then, the approximate eigenvalues of H(s) up to second order are given by 



(3-36) 



E^^s-l^-sE? + s 2 Ef\ 
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such that, the energy gap g(s) is given by 

g(s) = E X -E Q ^ 7 (* - l)(*i " *o) " «( =5 L T x " ^ L T + 

s 2 1 1 

7(^1) { S S - *;) E^=o ^0^) " S S - E^=o ^ 2 (^o)P/ 2 (^) } ' (3 " 37) 
Denoting the terms in the parentheses of the second term in (3-37) by A and the terms in the 

bracket of the last term by B, the energy gap is written as: 

s 2 

g( S )^ 1 (s-l)(x 1 -x )-sA+- —B. (3-38) 

7(s - 1) 

Now, in order to obtain the critical value of 7 in which the gap g(s) is minimum, we take the 
derivative of g(s) with respect to 7 so that = 0. Then, one can obtain 



7 - =rW(^' S#0 ' L (3 " 39) 



By substituting 7 cr4t in (3-38), we obtain 



g mm {s) = -s[A + 2JB{ Xl -x G )]. (3-40) 



Now, by using (3-26), the search time is evaluated as 

ds 



t 



"A 



es 2 (A + 2y/B( Xl - x )) 2 es(A + 2^B( Xl - x )) 2 
By substituting s — 1, the total search time at which the marked state \m) is found, is given 
by 

T = 1 (3-41) 

e(A + 2 v / J B(a; 1 - x )) 2 



As the above result indicates, in order to obtain the search time, one needs to evaluate the 
terms A and B defined below of (3-37) in terms of the corresponding graph polynomials Pi(x). 
On the other hand, A and B can be calculated easily via a systematic computer program for 
calculating the corresponding polynomials from the three term recursion relations (2-11). 
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4 Examples 

1. Dihedral group graph 

The dihedral group G = D 2n is generated by two generators a and b as follows: 

D 2n =< a,b : a n = 1, b 2 = 1, b~ l ab = a' 1 > . 

We consider the even n = 2m, where the odd n can be considered similarly. In this case, the 
conjugacy classes of the group, are given by 

C = {e}, Q = {a\ a- 1 }; i = 1, m - 1, C m = {a m }, C m+1 = {a 2j b; < j < m - 1} 

C m+2 = {a 2 ^ 1 6;0<j<m-l}, 

so we have m + 3 conjugacy classes. It is well known that the adjacency matrices defined as 
M = H g £d -ft(flO> are correspond to the underlying graph of the so-called group association 
scheme [22, 23] associated with the group D 2n , where R(g) is the regular representation of the 
element g of the group [24]. Therefore, we have 

A = I 2 ®I W Ai = I 2 ® (S i + S~ i ),i = 1,2,..., m- 1, 

A m = I 2 (g) S m , A m+1 = a x (g) (I n + S 2 + ... + S^), 
A m+2 = a x ®(S + S 3 + ... + S 2(m - 1 '>) 



where, S is the shift or circulant matrix defined as 

^01 
1 




V 



••. ••. ••. 1 

10 ••• 

and a x is the Pauli matrix. In order that we obtain a connected undirected graph, the adjacency 
matrices must be symmetrized, i.e., Ai = A\. To this end, we introduce the new adjacency 



Structured quantum search 



16 



matrix A' as A' = A rn+ i + A m+2 = cr x <g> (I + S + ... + S n x ) so that the Krylov bases are given 
by 

I0o) = {|e», 
1 



{|6> + |ba) + 16a"- 1 )}, 



102) = 



1 



and the adjacency matrix is represented as 



{| a ) + |^) + | a n-l)} 



A' 



n 





y/n \J n { n ~ 1) 

v yjn(n - 1) 



By replacing a = a,\ = a 2 = 0; f3i = y/n, f3 2 = \Jn(n — 1) in the three term recursion relation 
(2-13), the corresponding polynomials are obtained as 

Q (x) — 1, Q 1 (x)—x, Q 2 (x)=x 2 -n, 

Q 3 (x) = x(x 2 - n 2 ). 

So, the corresponding eigenvalues are given by x = —n, x± — and x 2 = n (the roots of 
Q 3 (x)). Then by replacing in (3-37), (3-40) and (3-41) respectively, and considering n >> 1, 
we obtain 

„ 2n-3 „ 4n-3 1 

A = — = 1, B — 



Icrit. 



2n 

s V4n-3 
1 - s 2n 2 

T ! 



4n 3 2n 2 



a • (s) - s( ^±V^E3^) 



2n 



An 2 



e(2n + Ay/E) 2 ' 
2. m-partite Graph 

An m-partite graph (i.e., a set of graph vertices decomposed into m disjoint sets such that no 
two graph vertices within the same set are adjacent) such that every pair of graph vertices in 
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the m sets are adjacent. Considering an m-partite graph which has n vertices in each of its 
disjoint sets, the corresponding adjacency matrix is given by 

A = K m ® J n , 

where, K m is the adjacency matrix of the complete graph with m vertices and J n is n x n all 
one matrix. Now by using the Lanczos iteration algorithm, the Krylov bases are obtained as 



|0i 



\<h 



1 



(\n + 1) + \n + 2) H h \n(m - 1))), 



y/(n - 1) 



(|2> + |3> + ... + |n» 



where, the adjacency matrix is reduced to the following 3x3 tridiagonal form 

/ n r~, n \ 



o 



o 



V 



\Jn(m- 1) 

n{m — 1) n{m — 2) ^n(n — \){m — 1) 

y/n{n- l)(m- 1) 



/ 



Now, as illustrated in the previous example in detail, by replacing a = a 2 = 0, a\ = n{m — 2) 



and Pi = yn(m — l),/3 2 = yn(m — l)(n — 1) in three term recursion relations (2-13), and 
evaluating the polynomials Pi(x) for i — 0, 1, 2 and Qatx), one can obtain 

x = —n, Xi = 0, x 2 = n{m — 1) 



and 

m(n — 1) — 1 „ 1 — n A — (m — l) 2 , 1 r m 2 n-2m + l, 

A = — , B = TV + \^T^T( TV' 

mn n n z m{m — 1) nm n z m z {m — 1) 

where for n ^> 1, m ^> 1, A and B are approximately as 

A = l, B = ^-. 

n 2 

Then, by using (3-39), (3-40) and (3-41), one can obtain 



1 „\/„3^ V 
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/ yjj — ^ 

9min(s) = S(1 + 2W— ), 



n 



and 



respectively. 



£(1 + 2^)2' 



3. Crown Graph 

A crown graph on 2n vertices is an undirected graph with two sets of vertices Ui and Vi and 
with an edge from Ui to Vj whenever i ^ j. The adjacency matrix of this graph is given by: 



A = K n ®o x 



where, K n is the adjacency matrix of the complete graph with n vertices and o x is the Pauli 
matrix. Then, the Krylov bases are given by 
I0o) = |1>, 



I0i) = 



|n + l) + |n + 2) + -.- + |2n-l)), 



|0 2 ) = ^=(|2) + |3) + --- + |n)), 
10s) = |2n). 



In the above bases, the adjacency matrix is represented as 

/ 



\ 



Vn^T 

^fn^l ra-2 

n-2 Vh~^l 

v 7 ^^ 

Again, by using (2-13) and (3-37) one can easily calculate the roots Xi and the quantities A 
and B as follows: 

x = -(n-i), x ± = -1, x 2 = 1, z 3 = (n-l), 



V 



/ 
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A 



n-2 



2n 



B = 




Now, by using (3-39), (3-40) and (3-41), one can obtain 



Icrit. — 



1 - S 



s 



n — 1 

n 2 



gmin(s) = S(l + 2 



(n-2)(ra- l) 2 
n 2 



1 



1 




en 



5 Conclusion 

Based on the local adiabatic evolution of the Hamiltonian and some techniques such as Krylov 
subspace projection methods, Lanczos iteration algorithm and spectral distribution methods, 
a new procedure for investigating quantum search in the state space of a graph was introduced, 
where an approximate analytical formula for calculating the corresponding minimum energy 
gap and the total search time was given. 
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